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ABSTRACT 

We present submillimetre (SCUBA) photometry of a sample of eight high-redshift 
(2.5 < z < 3.5) radio quasars from two redshift surveys: the TexOx-1000, or TOOT, 
survey; and the 7C Redshift Survey (7CRS). Unlike the powerful high-redshift radio 
sources observed previously in the submillimetre, these radio sources are typical of 
those dominating the radio luminosity density of the population. We detect just two 
of the T00T/7CRS targets at 850 /im, and one of these detections is probably due to 
synchrotron emission rather than dust. The population represented by the other six 
objects is detected in a statistical sense with their average 850 /im flux density im- 
plying that they are similar to low- redshift, far- infrared-luminous quasars undergoing 
at most moderate (^200 Mq yr~^) starbursts. By considering all the SCUBA data 
available for radio sources, we conclude that positive correlations between rest-frame 
far-infrared luminosity Lfir, 151-MHz luminosity L151 and redshift z, although likely 
to be present, are hard to interpret because of subtle selection and classification biases, 
small number statistics and uncertainties concerning synchrotron contamination and 
iT-correction. We argue that there is not yet any compelling evidence for significant 
differences in the submillimetre properties of radio-loud and radio-quiet quasars at 
high redshift. 

Key words: galaxies: active - galaxies: evolution - galaxies: formation - galaxies: jets 
- galaxies: luminosity function, mass function 



1 INTRODUCTION 

The advent of the Submillimetre Common-User Bolometer 
Array (SCUBA; Holland et al. 1999) camera on the James 
Clerk Maxwell Telescope ( JCMT) opened up the field of sub- 
millimetre cosmology. This, and other, instrument /telescope 
combinations now routinely detect dusty star-forming galax- 
ies at high redshifts, a process greatly aided by the way in 
which the steep submillimetre spectral index compensates 
for distance-dependent dimming, so that two objects of simi- 
lar far-infrared luminosity have similar observed submillime- 
tre flux density if both lie anywhere in the redshift range 
1 < z < 10 (Blain & Longair 1993). 

A major problem in understanding how submillimetre 
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surveys constrain the history of star formation in the Uni- 
verse is the difficulty of measuring redshifts for the sources of 
the submillimetre emission. This problem has not yet been 
solved, with only small (< 10 per cent) fractions of objects 
in submillimetre-flux-density-limited samples having reliable 
spectroscopic redshifts (e.g. Blain et al. 2002). 

An alternative strategy for studying the cosmic star- 
formation history is to make targeted submillimetre obser- 
vations of objects with known redshifts (e.g. Hughes, Dunlop 
& Rawlings 1997). Powerful radio sources are important tar- 
gets, not because the amount of star formation associated 
with them seems likely to be cosmologically significant, but 
because of two crucial properties: (i) they appear to be asso- 
ciated with galaxies, and also presumably black holes, with 
a mass distribution with a large mean and a small dispersion 
(e.g. Jarvis et al. 2001b); and (ii) studies of the physics of 
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their extended radio structures allows the calculation of the 
time since their jet-producing AGN were first triggered (e.g. 
Blundell & Rawlings 1999). The first property is potentially 
very useful for determining how physical properties derived 
from submillimotro observations (e.g. dust mass and star- 
formation rate) depend on galax;y, or black hole, mass. The 
second property can be used to constrain models in which 
starburst and AGN activity are physically related. 

The first large study of the submillimetre emission of 
high-rcdshift radio sources was made by Archibald ct al. 
(2001). They used SCUBA on the JCMT to make sensitive 
(to an 850 /xm noise level w 1 mjy rms) observations of 47 
radio galaxies covering the rcdshift range 1 < z < 4.5. The 
chief result of their study was a rise in the detection rates 
of radio sources from « 15 per cent at « ;$ 2.5 to ^75 per 
cent at higher redshifts. They interpreted their results as 
indicating that the submillimetre luminosity of radio galax- 
ies, and hence physical quantities like dust mass, increases 
systematically with redshift. 

One important limitation of the Archibald et al. (2001) 
study was that it targeted only radio galaxies and excluded 
radio quasars^. This was largely a practical decision, be- 
cause it was thought that the bright radio-emitting cores of 
quEisars might produce troublesome synchrotron contamina- 
tion. The decision could be defended on the basis of the sim- 
plest unified scheme (e.g. Antonucci 1993) in which no differ- 
ences are predicted in the joint distributions of isotropically- 
emitted quantities like synchrotron low-frequcncy-radio and 
dust submillimetre emission. A follow-up study of the sub- 
millimetre properties of powerful radio quasars at z ~ 1.5 
(Willott et al. 2002a) revealed a significant difference be- 
tween the submillimetre luminosities of radio-luminosity- 
matched quasars and galaxies - the quasars are brighter 
in the submillimetre by a factor > 2, accounting for syn- 
chrotron contamination. Willott et al. show that their result 
can still be accommodated within a unified scheme provided 
one allows for systematic changes in the properties of the 
obscuring material with time after the jet-triggering event 
(e.g. Blundell & Rawlings 1999; Hirst, Jackson & Rawlings 
2003) and/or correlations between quasar luminosity and 
the probability of viewing a given object as a quasar (e.g. 
as in the 'receding torus' model; Lawrence 1991, Simpson 
1998). 

A second important limitation of the Archibald et al. 
(2001) study was that it was confined only to the most radio- 
luminous sources: such objects contribute marginally to the 
luminosity density of the radio source population and could 
plausibly be entirely unrepresentative. The main aim of this 
paper is to establish the submillimetre properties of typical 
high-redshift radio sources. 

In Sec. 2 we describe how a sample of radio sources was 
selected for SCUBA follow-up, and how these observations 

t We will use the term radio quasar to refer to lobe-dominated 
steep-spectrum radio sources associated with quasar nuclei. These 
are the dominant population of radio-loud quasars in the low- 
frequency-selected samples considered in this paper, whereas 
compact flat-spectrum quasars can dominate in high-frequency- 
selected samples. Unified schemes (e.g. Antonucci 1993) invoke 
an optically-thick obscuring torus which hides the quasar nucleus 
when the jet axis makes a large angle with the line-of-sight, and 
such objects are termed radio galaxies. 



were made. The results of these observations are presented 
and discussed. In Sec. 3 we use these results to establish the 
submillimetre properties of typical high-z radio sources, to 
look for evidence of correlations between far-IR luminosity, 
radio luminosity and rcdshift and to compare and contrast 
these results with those obtained previously by Archibald 
et al. (2001) and WiUott et al. (2002a). Some concluding 
remarks are made in Sec. 4. 

We use J2000.0 positions throughout. The convention 
for all spectral indices, a, is that flux density 5*,^ oc ™, 
where v is the observing frequency. We assume through- 
out a low-density, A-dominated Universe in which Ho = 
70 km s"^Mpc"\ = 0.3 and nA = 0.7. 



2 SCUBA OBSERVATIONS OF 

HIGH-REDSHIFT TOOT AND 7CRS RADIO 
SOURCES 

2.1 Selection of SCUBA targets 

Our SCUBA targets were chosen from two of the suite of 
redshift surveys of low-frequency-selected radio sources out- 
lined by Rawhngs (2002): namely, (i) the TexOx-IOOO, or 
TOOT, Survey (Hill & Rawlings 2003); and (u) parts I & 

II of the 7C Redshift Survey (7CRst; Willott et al. 2002b). 
We aimed to focus on the specific region in the 151-MHz 
radio luminosity Lisi versus redshift z plane illustrated in 
Fig. 1. By studying objects in the region of intersection of 
the lightly-shaded horizontal and vertical bands of Fig. 1 we 
are sensitive to the one dex range in Lisi which contributes 
~ 50 per cent of the radio luminosity density across the 
range in cosmic time t from 1.8 Gyr, aX z = 3.5, to 2.6 Gyr 
at 2 = 2.5 (Willott et al. 2001b). 

We used the following selection criteria: for TOOT (Hill 
& RawUngs 2003), the TOOTOO and TOOT08 sub-regions 
with 0.1 < Si5i/Jy < 0.2, 2.5 < z < 3.5 and (for a prelim- 
inary version of the survey) _R^23.5, yielding five targets; 
for 7CRS, Sisi/Jy > 0.5 and z > 2.5, yielding four poten- 
tial targets of which one, 5C6.291, was excluded from our 
SCUBA target list because, with a 1.3-mm (IRAM) flux 
density of 65 ± 6 mJy (Steppe et al. 1995), its submillime- 
tre emission is clearly dominated by a l)riglit syiiclirotroii- 
emitting core. All objects are radio quasars: for the TOOT 
sub-sample this is probably because any radio galaxies meet- 
ing the Si5i and z criteria would be too faint optically to 
meet the R selection criterion; for 7CRS it is because there 
are no spectroscopically-confirmcd radio galaxies in the tar- 
get region, and none of the few objects without secure red- 
shifts are likely to be at suitably high redshift (Willott, 
Rawlings & Blundell 2001a). Basic data on the combined 
TOOT/7CRS sample can be found in Table 1. 

■t Parts I & II of the 7CRS were pursued in regions of the sky 
covered previously by the 5C radio surveys and, for consistency 
with previous papers, we use the 5C names throughout. 
^ We deem an object to be hyperluminous if, after making the 
corrections for synchrotron contamination detailed in Table 2, it 
has a far-infrared (30 — 1000 /im) luminosity Lpip^ > 10^'^ Lq; 
such objects have a rest-frame infrared (1 — 1000 /xm) luminosity 
LiR S> 10^^ Lq and, assuming a starburst heats the dust, a star- 
formation rate SFR> 1000 M© yr"! (see Fig. 3). The nine hyper- 
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Figure 1. The 151-MHz radio luminosity -L151 versus redshift 2 plane for various complete redshift surveys of radio sources, details 
of which can be found in the following references: the 3CRR, revised 3C, sample (triangles) — Laiiig, Riley & Longair (1983); the 6CE 
samijle (circles) — Rawlings, Bales & Lacy (2001); the 7CRS, 7C Redshift Survey (squares) - Willott et al. (2002b); five objects from the 
TOOT survey (stars) - Hill & Rawlings (2003). Objects classified as radio quasars, following the prescription of Willott et al. (1998), are 
shown as filled symbols; radio galaxies are shown as open symbols, and these have slightly larger size in the few cases where the redshifts 
are non-spectroscopic. The intersection of the lightly-shaded horizontal and vertical bands is the region of this plane [2.5 < z < 3.5 
and 26.75 < logj()(Li5i) < 27.75] central to the work of this paper. The darkly-shaded band shows the strip in radio luminosity 
considered in detail by Archibald et al. (2001). The dashed curves show the loci of objects with Oj-ad = 1 and 5i5i = 0.1 Jy (bottom 
curve), Sisi = 0.2 Jy (lower middle curve), Sisi = 0.5 Jy (upper middle curve) and Sibi = 2.0 Jy (top curve). Large rings represent 
hyperluminous infrared galaxies associated with powerful radio sources^, with the smaller rings showing all those objects observed, but 
not necessarily detected, with SCUBA (Archibald et al. 2001; Willott et al. 2002b; this paper); note that this includes objects, like the 
very-low-Li5i object 53W069 (labelled), not in the complete redshift surveys. 



2.2 SCUBA photometry and origin of the 
submillimetre emission 

The sample of T00T/7CRS radio sources was observed in 
photometry mode with the SCUBA bolometer array at the 
JCMT on the nights 2001 February 24-26, and 2001 July 30- 
31. Observations were made simultaneously at 850 fim and 
450 /Ltm. Our aim was to reach an rms sensitivity of ~ 1 mjy 
at 850 fim if 110 clear (> 5a) detection was achieved earlier. 
The data were reduced using the SURF package, following 
the same methods as Archibald et al. (2001) and Willott 
et al. (2002a). The resulting submillimetre flux densities are 



luminous objects are: 30181, 3C298, 3C318, 30432, 601045+3513 
(Willott et al. 2002a); 4C60.07, 4041.17, 801435-1-635 (Archibald 
et al. 2001); TOOT08.061 (this paper). 



given in Table 1. Note that the error bars on the photometry 
do not include systematic uncertainties in the absolute flux 
calibration which, from inter- and intra-night variations, we 
estimate to be w 15 and w 30 per cent at 850- and 450-/Ltm 
respectively. 

Just two of the eight objects are securely detected at 
850 fj,m, and only one of the eight objects at 450 jj,m. In 
Fig. 2 we show the spectral energy distributions (SEDs) 
of the two detected objects. From the comparison of these 
SEDs with other well-studied objects, it seems clear that the 
submillimetre emission of the one object detected at 450 fim 
(TOOT08_061) is dominated by dust emission. The submil- 
limetre emission from the other detected object (5C6.288) 
is probably dominated by synchrotron emission, although a 
signiflcant contribution from dust cannot be ruled out. 

We have coadded the 850- and 450-/im flux densi- 
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Table 1. Basic data for the sample observed with SCUBA from Rawlings, Hill & Willott (2004) for the TOOT objects and Willott et 
al. (1998) for the 7CRS objects. Following the classification scheme of Willott ct al. (1998), all objects are quasars, but all show evidence 
for obscuration. Wc have therefore sub-divided them according to our best guess as to the nature of the obscuration: Q(R) means the 
quasar nucleus is seen in lightly-reddened, transmitted light; Q(S) means the quasar nucleus is probably revealed by scattered, rather 
than transmitted, light — no objects are classified Q(N), meaning a naked quasar nucleus. The absolute B magnitudes arc not corrected 
for any extinction. Radio core luminosities have been estimated as close as possible to a rest-frame frequency of 5 GHz; a '-' denotes that 
existing radio data are of too poor resolution to allow isolation of any radio core. A f denotes that the projected linear size D is that 
between the core and an extended component on just one side of the core. The submillimetre photometry comes from the programme 
described in Sec. 2.2: detections at the > 2cr level are given in bold type; no corrections have been made for possible contamination by 
synchrotron radiation. 



ties of the other six targets obtaining an w 2.5(7 'statis- 
tical' detection at 850 //m given by the variance-weighted 
mean 1.22 ± 0.48 mJy at 850 /xm, and data consistent with 
pure noise at 450 (corresponding to a 2a upper limit 
~ 9 mJy). Repeating this exercise for just the four TOOT 
objects without clear submillimetre detections produces sim- 
ilar results (1.17 ± 0.58 mJy at 850 /um), although the sig- 
nificance of the statistical detection at 850 /um falls to « 2a. 
We conclude that more TOOT and 7CRS objects would bo 
securely detected at 850 /Ltm with only slightly deeper pho- 
tometric observations, although we cannot rule out a 'bi- 
modal' distribution in 850-/im flux density Ssbo such that 
some sources would only be detected with extremely deep 
exposures. 



3 THE SUBMILLIMETRE PROPERTIES OF 
TYPICAL HIGH-REDSHIFT RADIO 
SOURCES 

3.1 Spectral Energy Distributions 

We consider first the SED of TOOT08_061 (Fig. 2), the 
one T00T/7CRS object to be detected at both 450 and 
850 /xm. We used the 850/Ltm detection to fix the normali- 
sation of a template quasar SED in the following manner. 
The model SEDs arc the sum of three components: (i) a 
'starburst' component modelled by an isothermal optically- 
thin grey body of temperature 40 K and emissivity index 
/3 = 2 (see Priddey & McMahon 2001 for a justification 
of this choice of parameters); (ii) a 'Scyfcrt' component, 
dominant at 3 — 30 /im from Rowan- Robinson (1995); and 
(iii) a quasar component, dominant in the optical and near- 
IR, again from Rowan- Robinson (1995). The frequency- 
integrated luminosities of the three components are fixed 
in the ratio 0.3:0.3:1 so as to reflect the typical values de- 
rived by Rowan- Robinson (1995). This SED is consistent 
with the 450^^m data point given that it is associated with 
large random and systematic errors. Accounting for host 



galaxy reddening, the model SED also fits well to the opti- 
cal luminosity of TOOT08_061. Integrating under the SED, 
we estimate a bolometric luminosity ~ 7 x 10^^ Lq, and 
since the rest-frame infrared (1 — 1000 /im) luminosity Lir 
totals ~ 4 X 10^-^ Lq, TOOT08_061 is the one object in 
the T00T/7CRS sample with sufficient dust emission to 
classify the object as a hyperluminous infrared galaxy (e.g. 
Rowan- Robinson 2000). 

The correlation between submillimetre luminosity and 
bolometric luminosity for distant radio quEisars (Willott et 
al. 2002a) makes it no surprise that the brightest submil- 
limetre source in our T00T/7CRS sample is also the most 
optically-luminous quasar (sec Table 1). The scatter in this 
correlation means that it is equally unsurprising that the 
other optically-luminous quasar, 5C6.95, is detected only at 
the « 2(7 level with SCUBA. 

One guide wc have to the more typical properties of 
high-redshift radio sources comes from an 'average' TOOT 
SED, calculated from simple arithmetic mean values at 
each frequency, excluding data on TOOT08_061. This SED, 
which is clearly biased to a somewhat lower normalisation 
by the exclusion of TOOT08_061, is also plotted in Fig. 2. 
Synchrotron contamination of the Ssso ~ 1 nijy statisti- 
cal detection is a slight concern, but, as is illustrated in 
Fig. 2, extrapolation from radio to sub- mm wavelengths with 
Qrad = 0.7 [the value for the core of 5C6.288, and seemingly, 
from the discussion of Willott ct al. (2002a), a slightly flat- 
ter value than is typical for lobe-dominated radio sources], 
does not provide a dominant contribution at 850 /xm even 
if all the NVSS flux density arises from the core. Again fix- 
ing the normalisation of a standard quasar SED by the 5*850 
value (Fig. 2), wc find consistency with the 450/im upper 
limit and, by integration of the SED, an infrared luminosity 
Lir ~ 5 X 10^^ Lq, dominated presumably by emission from 
dust. 

Accounting for a small amount of dust extinction 
(amounting to Ay = 0.5 of Milky- Way-type dust as is suf- 
ficient to explain the red colours and narrow-Lya lines of 



© 0000 RAS, MNRAS 000, 000-000 



Submillimetre photometry of high-redshift radio quasars 5 



- 14 

13 ^ 

12 

• I — I 

m 

O 

S 11 ^ 
B 

^ 10 

9 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

^ 5C6.288 ■ 



TOOT08_061 
Average TOOT 



9 



7 / 

I I 



I I 



/. j' . I 



I I I 



10 



11 



12 



13 14 15 16 



Logio(Rest Frequency v / Hz 



Figure 2. The spectral energy distributions (SEDs) of objects in the TOOT/7CR,S sample. Data points are shown by the following 
symbols: TOOT08_061 (circles); variance-weighted mean of the other four high-rcdsliift TOOT quasars (squares, labelled 'Average 
TOOT'); 5C6.288 (stars). Error bars represent just the random errors. Open symbols represent data points corrected for dust extinction 
(amounting to Ay = 0.5 of Milky- Way-type dust) in the host galaxy. The dotted lines show extrapolations of the radio flux densities 
assuming a spectral index of 0.7. The dot-dashed lines show model SEDs, normalised to the 850 /im points, and based on the average 
quasar SED of Rowan-Robinson (1995) as detailed in Sec. 3.1. 



the TOOT quasars; Rawlings ct al. 2004), the model SED 
looks a little too bright at optical wavelengths but, given the 
crudeness both of the averaging process and the dust extinc- 
tion correction, the diflForence between model and data is not 
very significant. Thus, to the very limited extent to which we 
can currently define a typical SED for high-redshift TOOT 
sources, they look similar in shape to those of the low- 
redshift quasars studied by Rowan- Robinson (1995). The 
normalisation inferred for the average SED of the TOOT 
quasars is such that, in terms of infrared luminosity, typical 
TOOT quasars lie comfortably within the spread exhibited 
by low-redshift, optically-selected quasars. They are nat- 
urally classified alongside low-redshift quasars which have 
similar far-IR luminosities to ultraluminous infrared galax- 
ies; about one-third of the quasars in optically-selected low- 
redshift samples have these properties (Rowan-Robinson 
1995). 

There is of course a continuing debate as to whether 
star formation or an AGN dominates the heating of the 
dust responsible for the far-infrared emission from ultralu- 
minous infrared quasars. We can add little to this debate, 



noting simply that integration of the average TOOT SED 
of Fig. 2 gives a rest-frame far-infrared (30 — 1000 /jm) lu- 
minosity LpiR ~ 2 X 10^^ L0 which, in the absence of any 
quasar heating of the cooler dust, corresponds to a Icss-than- 
extreme star-formation rate SFR ~ 200 Mq yr^'^ (e.g. from 
equation [1] of Hughes ct al. 1997). Given the possibility of 
quasar heating, this estimate of SFR should probably be 
taken as an upper limit. 



3.2 Correlations with L151 and redshift 

We have calculated variance-weighted mean values of Ssso 
in various bins of 151-MHz radio luminosity L151 and red- 
shift z using the 66 radio sources with SCUBA data plotted 
with the smaller rings in Fig. 1. We eliminated from consid- 
eration the two sources whose 850/Ltm data are dominated 
by synchrotron emission (5C6.288 and 6C0902-I-34), and we 
made two assumptions concerning synchrotron contamina- 
tion of the other objects; (i) zero synchrotron contamination 
at 850/im for all objects; (ii) contamination at the levels de- 
tailed in the caption of Table 2. These assumptions should 
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bracket the uncertainties due to synchrotron contamination. 
Table 2 contains the results of this analysis, and far-infrared 
luminosities Lfir calculated from these results are plotted 
against redshift z in Fig. 3. 

3.2.1 Correlation with Lxf,x 

The baseline in Lisi opened up at high redshift by 
the TOOT/7CRS objects aUows us to look directly for 
a correlation between Lfir and Lisi at early cosmic 
epochs. The results of a straight comparison of the eight 
T00T/7CRS objects with the six more radio luminous 
but similar-rodshift objects (8C1039+68, MG1016+058, 
4C24.28, 4C28.58, 6C0902+34 and 6C1232+39) from the 
Archibald et al. (2001) study, excluding the two synchrotron- 
dominated objects (5C6.288 and 600902+34), are included 
in Table 2 and Fig. 3. Over this restricted redshift range 
(2.5 < z < 3.5) we see that, despite a one dex difference 
in average Lisi , there is no difference in average Sszo ■ How 
does this constrain any correlation between Lisi and Lfir? 

First, we neglect synchrotron contamination and con- 
sider the effects of the submillimetre difference between 
radio quasars and radio galaxies discovered by Willott et 
al. (2002a). The average value of iSsso for the « ~ 3 
T00T/70RS objects is likely to be biased high by the ex- 
clusion of radio galaxies (sec Sec. 2.1), and the average Ssso 
for the highcr-Li5i z ~ 3 objects biased low by the exclusion 
of radio quasars (Archibald et al. 2001). We can estimate the 
relevant bias factors as follows. As the higher-Lisi sources 
have similar values of Lisi to the radio quasars studied by 
Willott et al. (2002a), it is reasonable to assume that they 
have a similar quasar fraction (« 0.4; Willott et al. 2000) 
and a similar submillimetre difference (radio quasars being, 
after correction for synchrotron contamination, a factor ~ 3 
brighter than the radio galaxies at 850 /im), so that the 
average Ssao would be expected to be a factor ~ 2 higher 
in a genuinely complete sample. For a genuinely complete 
TOOT/7CRS sample, the average Sszo is likely to be lower, 
but by a smaller factor because the quasax fraction is likely 
to be similar (Willott et al. 2000) and there will a lower con- 
tribution to the average by the excluded objects. Overall, 
therefore, for z ~ 3 complete samples, the high-Lisi point 
in Fig. 3 is expected to be higher than the low-Lisi point 
by a factor ~ 2. This is a similar separation to that seen 
between the high- and low-Lisi points at 2 ~ 1 and 2 ~ 2 
in Fig. 3, and these do not need large corrections for ob- 
ject exclusion because the Willott et al. (2002a) programme 
ensured that the SOUBA-targeted radio quasars and radio 
galaxies ended up in roughly the same ratio as the quasar 
fraction. 

We conclude that, although significant only at the ~ 2cr 
level (see also Sec. 3.2.2), there is evidence for a Lfir-Lisi 
correlation in Fig. 3. An increase in Lfir by a factor ~ 2 as 
Li5i increases by a factor ~ 10 are consistent with a model 
in which Lfir oc Lfgj, provided p^O.5. 

If we adopt the synchrotron-corrected values from Ta- 
ble 2, then a somewhat different picture emerges (Fig 3). 
Most noticeably, there now exists a large (factor ~ 10) dif- 
ference between the high- and low-Lir,i objects at z ^ 1. 
Synchrotron corrections have most effect at low redshift be- 
cause 850/ttm corresponds to a lower rest-frame frequency, 
and most of the SOUBA-targeted galaxies are from the 



radio-bright 30 survey". This might be interpreted as evi- 
dence that we have, if anything, underestimated the strength 
and slope of the Lfir -L151 correlation, or alternatively that 
the (uncertain) synchrotron corrections have been too se- 
vere. However, it is also not completely implausible that the 
synchrotron corrections are not severe enough [although see 
Willott et al. (2002a) for arguments to the contrary]. If this 
were the case, most 850-/im detections at « ~ 1 could be 
dominated by synchrotron contamination and the difference 
in dust submillimetre luminosity between radio quasars and 
radio galaxies inferred by Willott et al. (2002a) might have 
to be reduced. The correct interpretation of Fig. 3 might 
then be that there is little evidence for a Lfir-Lisi correla- 
tion at any redshift, but there is a dramatic drop in Lfir at 
2 ~ 1 (see Fig 3); we return to this point in Sec. 3.2.2. 

3.2.2 Correlation with redshift 

The similarity between high-a T00T/7CRS quasars and 
some low-redshift optically-selected quasars, discussed in 
Sec. 3.1, suggests we investigate whether there is any di- 
rect evidence for the cosmic evolution in the rest-frame far- 
infrared luminosities Lfir of radio sources of 'typical' radio 
luminosity, i.e. across the lightly-shaded horizontal strip in 
Fig. 1. The results of this investigation are included in Ta- 
ble 2 and Fig. 3. 

Again, we start our discussion by adopting the assump- 
tion of negligible synchrotron contamination. The evidence 
for any systematic rise in Lfir with redshift is then weak 
(< 2a) and, as was the case for the hints of a correlation 
between Lfir and L151 (Sec. 3.2.1), driven by the low aver- 
age Lfir of the I0W-L151 , 2 ~ 1 bin. Mindful of the many 
bieises discussed in Sec. 3.2.1, we conclude that, neglecting 
corrections for synchrotron contamination, there is no firm 
evidence for a cosmic rise in Lfir for typical radio luminos- 
ity sources. As we have already noted in Sec. 3.2.1, adopting 
the synchrotron-corrected values from Table 2 means that 
there is much stronger evidence for cosmic evolution in the 
sense of a dramatic drop in Lfir at z ~ 1. 

Evidence for any cosmic rise in Lfir for the high-Lisi 
sources (Fig. 3), with our without synchrotron correction, is 
very weak, driven, this time, by a high average Lfir in the 
2 ~ 4 bin, a redshift range not yet probed by SCUBA obser- 
vations of I0W-L151, i.e. TOOT, sources. There arc just three 
objects (60*0032-1-412, 801435-1-635, 60*0140-1-326) in this 
bin, all from the Archibald et al. (2001) SOUBA survey, so 
small number statistics and selection biases are extremely 
important. These objects are from 'filtered', i.e. steep- 
radio-spectrum-sclcctcd, rather than complete-flux-density- 
limited, samples (80, Lacy et al. 1994; 60*, Jarvis et al. 
2001a). For such samples, 60* is the only one in which it has 
been demonstrated that potentially high-z objects have not 
been missed because of low optical emission line strengths 
(Jarvis et al. 2001a). Given that scaling relations between lu- 
minosities in different wavebands are normally positive, the 



II Note that the method of synchrotron correction adopted by 
Archibald et al. (2001) tends to set small positive values of Sgso 
of low-^:, typically 3C, sources to zero, whilst leaving negative 
values unchanged. This may result in 'corrected' values which are 
biased low. 



© 0000 RAS, MNRAS 000, 000-000 



Submillimetre photometry of high-redshift radio quasars 7 



Sample 


Range in 






O860 at z 


— ± 


O8B0 at z ■ 




O860 at z - 


_ Q 
- O 


O860 at z - 


- A 


All in Fig. 1 


26.75 < Lisi 


< 


27.75 


n 


0.85 ± 0.43 


(5) 


1.87 ±0.35 


(9) 


2.40 ± 0.44 


(7) 


- 












y 


0.12 ±0.43 


(5) 


1.69 ±0.35 


(9) 


2.40 ± 0.44 


(7) 






All in Fig. 1 


27.90 < Li5i 


< 


28.50 


n 


2.52 ± 0.36 


(8) 


2.97 ±0.26 


(16) 


2.63 ±0.52 


(3) 


5.84 ±0.59 


(3) 










y 


1.55 ±0.36 


(8) 


1.82 ±0.26 


(16) 


2.47 ±0.52 


(3) 


5.77 ±0.59 


(3) 


Archibald ct al. 


27.90 < Lisi 


< 


28.50 


n 


1.80 ±0.41 


(6) 


2.10 ±0.30 


(11) 


2.63 ±0.52 


(3) 


5.84 ±0.59 


(3) 










y 


0.98 ±0.41 


(6) 


0.97 ± 0.30 


(11) 


2.47 ±0.52 


(3) 


5.77 ±0.59 


(3) 


Archibald et al. 


AH Li5i 






n 


0.96 ±0.27 


(13) 


1.60 ±0.22 


(21) 


2.82 ± 0.41 


(5) 


7.38 ± 0.38 


(7) 










y 


0.20 ±0.27 


(13) 


0.89 ± 0.22 


(21) 


2.70 ± 0.41 


(5) 


7.32 ± 0.38 


(7) 



Table 2. Variance-weighted mean values of 850-/um flux density 58bo for samples of radio sources covering various ranges of 151-MHz 
luminosity -Lisi and redshift z. These are calculated for both all the SCUBA-observed radio sources shown in Fig. 1 of this paper, and 

for just those objects in the Archibald ct al. (2001) study. Synchrotron-dominated objects (5C6.288 and 6C0902+34) have been excluded 
from the analysis, and values calculated both with (SC 'y') and without (SC 'n') synchrotron correction arc given. The synchrotron- 
corrected values are taken from Archibald ct al. excepting for the objects from Willott ct al. (2002a) and this paper. For the Willott 
et al. objects, we assumed a radio-to-submillimctrc spectral index a = 0.7, together with other constraints discussed by Willott ct al., 
to estimate the following synchrotron contaminations, expressed as fractions of the observed Ss5o- 3C181 (0.06); 3C191 (0.30); 3C205 
(0.52); 3C268.4 (0.58); 3C270.1 (0.40); 3C280.1 (0.35); 3C298 (0.38); 3C318 (0.00); 3C432 (0.05); 6C0955±3844 (0.00); 6C1045±3513 
(0.05). For the TOOT and 70 objects we estimate zero synchrotron contamination because such low-5i6i objects can only contribute 
significantly at 850 /um if their radio flux is concentrated in a compact core (see Fig. 2). The redshift bins used were 0.5 — 1.5, 1.5 — 2.5, 
2.5 — 3.5 and 3.5 — 4.5. The number of objects involved in calculating each average are given in brackets. 



use of samples in which low-emission-line-strength objects 
have been missed would have the effect of biasing the aver- 
age 5850 value high. Considering just the two 6C* objects, 
wc see in Fig. 3 no evidence that Lpir is any higher at z ~ 4 
than at lower redshift. A rise in average Lpir only emerges 
if one adds in the hyper luminous object 8C1435+635. We 
caution against putting much weight on either this rise, and 
its lower formal error bar, because the parent (8C) sample 
docs not have a well-defined redshift completeness. It should 
also be noticed that any 2; ~ 4 radio quasars are likely to 
have been missed by all filtered samples, including 6C*, be- 
cause their prominent radio cores can mask steep-spectrum 
extended emission. The inclusion of such objects would tend 
to increase the average value of Ssso (Willott et al. 2002a). 

Adopting the assumption of negligible synchrotron con- 
tamination, the weakness of any correlations with redshift in 
Fig. 3 means we should re-examine the seemingly strong ev- 
idence for a cosmic rise in Lfir from the study of Archibald 
et al. (2001). Binned data from their study are presented in 
Table 2 and plotted in Fig. 4. Considering just those objects 
in the high-Li5i strip of Fig. 1 we see that, as in Fig. 3, 
evidence for a cosmic rise comes from just the a ~ 4 bin, 
and is subject to the same worries concerning small-number 
statistics and selection biases. It is only by averaging all the 
Archibald et al. data, regardless of L151, that a seemingly 
monotonic rise in Lfir with z emerges (Fig. 4). There are 
several potential problems with interpreting this correctly. 
First, there is a strong and systematic increase in the av- 
erage L151 of objects from the lowest-redshift bin (mean 
Lisi = 27.5) to the highcst-rcdshift bin (mean L151 — 28.6): 
this could account for nearly all of the 'cosmic' rise in Lfir 
given a correlation of the form Lfir oc Lf gj , with as 
discussed in Sec. 3.2.1. Second, the average value of the z ^ 4 
bin is, despite the low formal error bar, calculated from ob- 
jects drawn from various filtered surveys, and therefore sub- 
ject to the possible selection biases discussed in Sec. 3.2.1. 
Third, the reality and magnitude of any cosmic rise in Lfir 
is strongly dependent on the assumed dust template. Our 



choice of dust template (dust temperature T = 40 K, emis- 
sivity index /3 = 2.0) is motivated by fits to the SEDs of 
distant quasars by Priddey & McMahon (2001). For this 
template, objects of a given Lfir have higher Sgr,o at z ~ 4 
than at ~ 1, whereas, as illustrated in Figs. 3 and 4, the 
Archibald et al. choice of dust template (T = 40 K, /3 = 1.5) 
means that at a given Lfir, 5850 is virtually independent of 
redshift. Thus, for a given cosmology and a given dataset, 
our choice of dust template will always produce less cos- 
mic evolution than the template adopted by Archibald et 
al. (2001). 

It is also worth re-emphasising that the Archibald et al. 
(2001) study excluded radio quasars. It could be argued that 
this reduces, or even eliminates, the need for the complicated 
corrections for biases discussed in the context of Fig. 3 in 
Sec. 3.2.1. However, radio galaxies will only show the same 
relationship between Lfir and redshift as the total (radio 
galaxy plus radio quasar) population if the probability of a 
radio source being classified as radio galaxy is independent of 
redshift. This is an as yet unproven assumption, and one that 
does not sit easily with the results of Hirst et al. (2003) and 
Willott et al. (2002a) which suggest that both the quasar 
fraction and submillimetre luminosity are functions of radio 
source age which, in turn, is likely to correlate negatively 
with redshift (e.g. Blundell & Rawlings 1999). 



4 CONCLUDING REMARKS 

We have seen in Sec. 3.2.2 that evidence for evolution of 
the dust properties of radio sources with redshift emerges 
only with the adoption of the synchrotron-corrected values 
for 5850, and is confused by the influence of the correlation 
between Lfir and L151 (Sec. 3.2.1). Since these corrections 
and this correlation arc very uncertain, we conclude that, for 
radio sources, direct observational evidence that Lfir rises 
systematically with redshift is not yet compelling. Although 
it would be rather surprising if an increase in the dust and 



© 0000 RAS, MNRAS 000, 000-000 



8 Rawlings et al. 



All SCUBA — observed radio sources 



1.00 ^ 



CO 



.10 



0.01 



27.90 ^ Logio(L,5,) ^ 28.50 
26.75 ^ Log,„(Z„„) ^ 27.75 



r = 


40 


K, 


/5 = 


2.0 - 


T = 


50 


K, 




2.0 


T = 


40 


K, 


/? = 


1.5 - 



-- 1000 



■CD" 

o 




J I L_ 



10 







2 3 

Redshift z 



4 



5 



Figure 3. The variance-weighted mean far-infrared (30 — 1000 ^in) luminosity LpiR of radio sources in selected ranges of 151-MHz 
luminosity Lisi versus redshift z (in bins 0.5 — 1.5, 1.5 — 2.5, 2.5 — 3.5 and 3.5 — 4.5). All objects identified as observed with SCUBA 
in Fig. 1, barring three synchrotron-dominated objects, were used to calculate these averages (see Table 2). We have calculated LpiR 
using an isothermal optically-thin grey body spectrum characterised by a temperature T = 40 A" and emissivity index i3 = 2; the locus 
of an 5850 = 1 mjy source with this spectrum is shown by the solid line, with the effects of varying T and j3 being shown by the dotted 
line (T = 50, /? = 2) and the dashed line (T = 40, /3 = 1.5). Lfir is converted to star formation rate SFR using eqn. (1) of Hughes 
et al. (1997), adopting e = 1. The large open squares represent radio sources with 26.75 < log]^g(Li5i) < 27.75. The large open circles 
represent radio sources with 27.9 < logiQ(Li5i) < 28.5, i.e. the high-luminosity strip considered by Archibald et al. (2001); two different 
values are plotted at z = A, the lower point representing just the two 6C* objects, the upper point representing all three objects in this 
bin. The smaller symbols, with error bars not plotted, show the mean values after correction for synchrotron contamination (see Table 2). 
Note that this plot is affected by the various potential systematic biases discussed in Sec. 3. 



gas masses of active galaxies with redshift were not an im- 
portant part of the overall picture (Archibald et al. 2001), 
the results of Sec. 3.2 arc now in reasonable accord with 
the results obtained on radio-quiet quasars. Priddey et al. 
(2003) noted that variations in the submillimetre detectabil- 
ity of optically-selected quasars with redshift are consistent 
with the K-correction of a template spectrum which does 
not evolve in normalisation with redshift. Although, sub- 
millimetre studies of radio-quiet quasars do not suffer from 
the problems of synchrotron contamination, it is interesting 
to note that, for radio-quiet quasars, searches for system- 
atic changes in dust emission with redshift are confused by 
possible correlations between ipiR and optical luminosity 
(Willott, Rawlings & Grimes 2003; Priddey et al. 2003). 

In Sec. 3.2 we have seen that hyperluminous infrared ob- 
jects contribute substantially to the average values of Sgso 
used in the correlation studies. Looking at the location in 



Fig. 5 of the radio-loud objects with associated hyperlu- 
minous infrared galaxies hints at complicated physical pro- 
cesses behind any correlations. As has been shown previously 
by Willott et al. (2002a), and was conjectured by Blundell 
& Rawlings (1999), there seems to be a clear association 
of hyperluminous infrared galaxies with the radio sources 
of smallest projected linear size D. The underlying physics 
here (see caption of Fig. 5 for details) may be that a short 
time r ^ 10^ yr since the jet-triggering event is a necessary, 
if not sufficient, condition for far-infrared hyper luminosity 
as at later times dust may have been destroyed or dispersed 
by the expanding radio source (Dc Young 1998). With two 
important exceptions, the hyperluminous objects appear to 
be concentrated amongst the most radio-luminous sources. 
The underlying physics here may be that extreme radio lu- 
minosity is strongly linked to the presence of an extreme- 
luminosity quasar nucleus which in turn promotes, perhaps 
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Figure 4. The variance-weighted mean far-infrared (30 — 1000 fim) luminosity LpiR of radio sources in selected ranges of 151-MHz 
luminosity L151 versus redshift z (in bins 0.5 — 1.5, 1.5 ~ 2.5, 2.5 — 3.5 and 3.5 4.5). Only objects included in the Archibald et al. 
(2001) study were used to calculate these averages (sec Table 2). The axes and lines have the same meaning as in Fig. 3. The large stars 
represent just the Archibald ct al. (2001) objects in the 27.9 < log^Q(Li5i) < 28.5 strip (two values arc plotted at z = 4, the lower point 
representing just the two 6C* objects, the upper point representing all three objects); the large triangles represent all the Archibald et 
al. (2001) objects, irrespective of radio luminosity. The smaller symbols, with error bars not plotted, show the values after correction for 
synchrotron contamination (see Table 2). Note that this plot is affected by the various potential systematic biases discussed in Sec. 3. 



by direct heating of dust, infrared hyperluminosity. The two 
hyperluminous objects at relatively low Lisi (6C1045+3513 
and TOOT08_061) both show, unusually for radio-loud ob- 
jects, evidence for broad absorption lines (Willott et al. 

2002a; Rawlings ct al., 2004). The underlying physics here 
may be that young bursts of less-than-extreme radio lumi- 
nosity, probably associated with the early evolution of less 
powerful jets (e.g. Kuncic 1999), as well as BAL features 
may be indicative of processes triggered by galaxy mergers 
(e.g. Canalizo & Stockton 2001). 
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